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ABSTRACT 

We extend previous work on sum rules for the invariant amplitudes of pion 
Compton scattering by deriving a complete lowest order perturbative spectral func- 
tion - and its leading non perturbative power corrections - for a specific combination 
of the two helicities (Hi + H 2 ) of this process. Using some properties of a modified 
version of the Borel transform, we develop a method of calculation of the gluonic 
corrections which can be easily extended to other similar reactions, such as proton 
Compton scattering. A preliminary comparison of the new sum rule with the pion 
form factor sum rule is made. 
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1 Introduction 



Among the semi-phenomenological methods developed in the past in the descrip- 
tion of the resonant region of QCD, sum rules are the closest to a fundamental 
description of the strong interactions and they have been formulated within a cal- 
culational scheme which is very similar to perturbative QCD. In the sum rules 
approach a "master equation" relates the (timelike) region of the external invari- 
ants in the correlation functions, the region where the resonances are located, 
to the euclidean region where information from the operator product expansion is 
parameterized by local vacuum condensates (power corrections). In practical appli- 
cations, only the phenomenological input from the lowest dimensional condensates 
(quark and gluon condensates) is required while the resonance region is described 
by a model spectral density of states which contains non perturbative information 
on the lowest resonance, although mixed with the continuum contribution. 

A particular class of sum rules are Borel sum rules. They are obtained by 
acting with Borel transforms on the master equation, in order to enhance the 
resonant contribution respect to the contribution from the higher states, while 
a superconvergence condition on the same spectral density allows to remove the 
continuum contribution beyond a given "duality interval" (sq). The net result of 
this procedure is a "sum rule" which relates properties of the lowest resonance to 
an operator product expansion (ope) series which includes both perturbative and 
non perturbative terms. This way of obtaining information about the resonances, 
which works very well for both 2 and 3 point functions, has largely contributed 
to our understanding of QCD in the intermediate energy region |], [§ [J. In the 
extension of the method from 3 to 4 point functions various new difficulties are 
encountered, such as 1) a more complex dependence of the dispersion relations on 
both the invariants s and t, and 2) the question related to the validity of the ope. 
As discussed in Ref . || , this last fact constraints the size of the kinematical s and 
t region in which the ope is applicable, demanding to consider by this method scat- 
tering at fixed angle {—t/s fixed). Once these kinematical conditions are satisfied, 
the new information available from sum rules for 4-point functions, compared to 
3 point functions, comes from the more general s and t dependence of the residui 
at the double poles of the lowest resonance (the invariant amplitudes), which can 
then be related to a perturbative spectral function and to its power corrections. 
Therefore, compared to form factors 0, |!| , in the case of 4-point functions a true 
angular dependence of the sum rule comes into play and a parallel study of the 
result obtained by this method and by the factorization picture, for Compton scat- 
tering, becomes possible. This last aspect has been preliminarly discussed in Ref. 
H], where predictions from the asymptotic local duality sum rule of Ref. || and a 
modified factorization formula which resums Sudakov effects have been compared. 

In the first part of this work (sections 2 and 3) we give a complete - i.e. non 
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asymptotic - derivation of the local duality sum rule for pion Compton scattering, 
extending the results obtained in Ref. |J and we show that the local duality sum 
rule can be organized in a rather compact form in the Breit frame of the pion, where 
the spectral densities become polynomial. For this purpose we use the cutting 
rules for the propagators, as discussed in |J. In sections 4 we use a modified Borel 
transform, first introduced in Ref. H, to evaluate the gluonic power corrections 
to the spectral function. An original scheme of calculation of such contributions 
is discussed. A more technical discussion of its implementation in the calculation 
of a specific coefficient of the expansion is briefly described in appendix B. In 
section 5 we discuss the evaluation of the power corrections proportional to the 
quark condensates and give the complete expression of the spectral funtion which 
appears in the sum rule for the sum [H\ + H 2 ) of the two helicity amplitudes of 
Compton scattering. In sections 6 the new Borel sum rule is derived and compared 
to the corresponding one for the form factor case, while our conclusions are in 
section 7. Our convention can be found in appendix A; section B briefly illustrates 
the evaluation of a specific gluonic correction. 

2 Overview 

We start our discussion by first quoting the main results of Ref. [|[ . As a starting 
point we consider the following four-point amplitude 

T afMuX (pl,pl, s,t) = i J d 4 xd 4 yd 4 zexp(-ip 1 ■ x + ip 2 ■ y - iqi ■ z) 

x(0\T(7i.(y)Jp(z)J v {0)7i{(x))\0) , (1) 

where 

^ = 2/3147^-1/3^7^, ri a = u^ a d (2) 

are the electromagnetic and axial currents respectively of up and down quarks. 
The choice of this specific correlator is motivated by the particular process that we 
have in mind, namely Compton scattering of a pion. The pions lines of momenta 
Pi and P2 are assumed to be off-shell while the photons lines (qi, q 2 ) are on-shell. 

The perturbative expansion of the correlator in eq. (1) is represented by the 
diagrams in Figs. 1-3. Specifically, Fig. 1 represents the perturbative contribution 
to the expansion, while Figs. 2 and 3 describe the power corrections which are 
proportional to the condensates of quarks (Figs. 2) and gluons (Figs. 3). They 
parametrize the non perturbative structure of the QCD vacuum. As discussed in 
Ref. pi, in the case of pion Compton scattering, when the two external photons 
are physical and all the Mandelstam invariants are moderately large (see also the 
discussion in Ref. ||), the contribution to the leading scalar spectral function 
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comes from the diagram in Fig. 5a. This diagram differs from the corresponding 
one in the form factor case (Fig. 5b) for having one additional internal line (the top 
line connecting the two photons, see Fig. 5a) off-shell. Notice that such contribu- 
tion is not related to any gluonic exchange, differently from the usual perturbative 
expansion based on factorization theory Q where the exchange of one gluon is 
supposed to be dominant. The main argument of Ref. || is that it is possible to 
write down a sum rule for Compton scattering in terms of a perturbative spectral 
function calculated (to lowest order) from Fig. 5a and from its power corrections. 
On the other hand, the resonance region is modeled, as usual, in terms of a single 
resonance (the double pole at m 2 , the pion mass squared) and a continuum contri- 
bution starting at a virtuality sq. Information on the invariant amplitudes of pion 
Compton scattering can be obtained from the residue at the pion double pole of 
the correlator in (1). 

In fact, the invariant amplitudes of pion Compton scattering can be isolated 
from the matrix element 

M uX = ij d 4 ye- i ^( P2 \T(J u (y)J x (0)) \ Pl ). (3) 

In Ref. we have expanded (3) in terms of Hi and H 2 defined by the relation 

AfV = Hi(s, t)e (1)A e (1) ^ + H 2 (s, t)e (2)A e (2)A \ (4) 

where and are polarization vectors, which are defined below. The residue 
at the pion pole is isolated from the exact spectral density of (1) as 

A^vx* = flvxuVio (2tt) 2 5(p 2 i - ml)5{pl - ml) 

xM uX (pi,p 2 ,qi) + A cont +h.r. (5) 

where M uX is given by eq. (4), A cont is the double discontinuity of the multiparticle 
states of the correlator in eq. (1), and h.r. denotes the contribution due to the 
higher resonances. In the derivation of the sum rule it is assumed that other 
single particle states, corresponding to the double poles located at momenta p\,p 2 
equal to the masses of higher resonances, give negligible contributions. Therefore 
ApuXcr contains the most valuable non-perturbative information on the invariant 
amplitudes Hi and H 2 at the pion pole, although mixed with the contribution of 
the higher states. While the contributions from A cont - for virtualities bigger than 
a threshold value so - cancel from both sides of the sum rule (see eqs. (10) and (11) 
below), a suitable way to remove the contribution of single particle states from (5), 
still keeping the pionic contribution (Hi and H 2 ), is to use a projection. In Ref. 
H this has been obtained using two projectors, defined in terms of the vectors 
e A' '> e i an d n \i the last one satisfying the requirements 

n ■ (qi) = n - (pi- p 2 ) = n 2 = 0, (6) 
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and possibly choosen as 

n»= (eW±ieV>Y . (7) 

Improvement of the sum rule is obtained, as usual, by acting on it with double 
Borel transforms |], ||, |j . In the case of pion Compton scattering the two projected 
spectral densities are given by the expressions 

A^ 2 \p 2 ,s,t) = f 2 w n-p 1 n-p 2 (2n) 2 5( P 2 1 -m 2 w )5(p 2 2 -ml) 

xiH^ + H^t)) +Ai 12)co " 4 , (8) 

and 

A^(p 2 ,s,t) = fln-p l n-p 2 {2^) 2 5{p\-ml)5{pl-ml) 

xflifot) +Ai 1)con *. (9) 

In the following we shall assume m n to be zero. Then the full sum rule or "master 
equation" for the invariant amplitudes of pion Compton scattering, for instance 
for Hi + H 2 , can be written down as 

f 2 n- Pl n-p 2 (ifi(s,t) + H 2 (s,t)) 

= fda! [ S °ds 2 A^ ert ( Sl ,s 2 ,s,t)e-^ M ?-^ 
Jo Jo 

(l _ e -(*si)/M^ ^ _ e -(A2- S2 )/M|^ + pQwer corr _ ; (1Q) 

where A 2 = (s+t) /4 is the radius of the contour in the spectral representation of the 
correlator (1) (see Fig. 4) and the remaining terms in (10) are the power corrections. 
In the local duality limit of eq. (10) the power corrections are neglected. Keeping 
M 2 and X/M 2 to be large, the sum rule for Hi + H 2 takes the approximate form 

U 2 n-pm-p 2 (Hi(s,t) + H 2 (s,t)))= r dsi f S ° ds 2 A^ rt (si,s 2 ,s,t), (11) 

Jo Jo 

in which the exponentially suppressed terms present in eq. (10) have been re- 
moved. A similar sum rule can be written down in terms of only Hi and A( 1 ) per *. 
The asymptotic expansion at fixed angle of the sum rule for the two perturbative 
spectral densities A^ pert and A( 12 ) pert , calculated from the ope-side of the sum rule, 
has been given in jj. In the next section we are going to discuss their complete 
evaluation to lowest order in a s . 



3 Perturbative spectral functions 

In this section we present the complete evaluation of the lowest order perturbative 
spectral functions which appears in the sum rule for the sum of the two helicity 
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amplitudes eq. fllPf) and for the single invariant amplitude Hi(si, S2, s, t). In a more 
complete analysis of the sum rule for pion Compton scattering, an asymptotic 
expansion of the lowest order spectral functions is no longer sufficient. In fact, 
while the dominant 0(1/Q 4 ) and the next to leading contributions can be obtained 
in a more simple way, by a simple power counting of the various integrals, the 
determination of the 0(1/Q 8 ) contributions requires a full scale calculations of 
all the terms from Fig. la. This turns out to be a rather lenghty task. The 
asymptotic expansion - arrested to 0(1/Q 4 ) - of the sum rule for H1 + H2 has been 
already given in Ref. ||. In that work it has been shown that, at fixed angle, 
Compton scattering has a strong similarity to the form factor case. This crucial 
observation has opened the road to the extension of the sum rule programme to 
4-point functions for Compton processes, while providing a calculational scheme 
which is a direct generalization of those known in the form factor case || . 

The leading spectral functions A^ pert and A( 12 ) per * can be calculated in terms 
of 3-cut integrals of the form (see Fig. 5a) 

ime, k . P1 , ...)] = id'km k- P1 , ...) ^)M( Pl -*)-)Mfe-*)') 

J (v 1 -k + 01 Y 



I'[f(k 2 , k ■ Pl , ...)]= / d 4 kf(k 2 , k ■ Pl , ...)M* 2 )M(Pi - kf)5 + ((p 2 - A;) 2 ), (13) 
The analogy to the form factor approach (see Fig. 5b) is clearly due to the fact 



that (|12[) approximates the total discontinuity (to lowest order in a s ) along the 
positive si, S2 virtualities only by a "3-particle cut", while one propagator (the top 
line of Fig. 5a) is kept off shell. The dominance of this discontinuity compared to 
other contributions can be easily proved. Other cuts are 1) either sub leading or 2) 
they are relevant for forward scattering (not considered in this work) or 3) they 
are not in the physical region. The evaluation of these integrals is not obvious and 
can be handled in a suitable reference frame. It simplifies considerably in the pion 
"brick- wall" frame in which the momenta of the two pions are parameterized as 

Pi = Qv + ^v, p 2 = -^v + Qv, (14) 

where v and v are light-cone momenta, such that v 2 = v 2 = 0, v ■ v — 0. The 
momentum transfer t can be expressed in terms of Q 2 , which acts as a large "pa- 
rameter" in the scattering process, by the relation 

t=(P2-Pi) 2 = s 1 + S2-2Q 2 - S 0. (15) 



Q 2 is given by 



Q 2 = ^{si + s 2 — t + 5) = -(s + u + 5) , (16) 
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with 

4Q 4 - s lS2 
2Q 



Sl + s 2 - ty - A Sl s 2 = ^ ^ 2 . (17) 



In the frame specified above, all the transversal momentum is carried by the two 
photons. Notice that t = —2Q 2 for si = s 2 = 0. In our approach the two photons 
are on shell (qf,q 2 = 0) and carry physical polarizations. It should be noticed 
that the matrix element in (3) admits various possible extensions (in terms of 
polarization vectors) when we move off shell the two pions. This extension must 
be done in such a way that the requirements of physical polarizations for the two 
photons are still satisfied. A suitable choice of the polarization vectors drastically 
simplifies the calculations. The two polarization vectors that we are going to use 
in our calculations are slightly different from those given in [[J. In our case they 
are defined as 

e (l)A = iV e (2)A _ £ h 



where 
and 



N x = e x ^P^r u R p P x = v lP x + v 2V \ + |# A , (19) 

v\=V\'V2-s 2 u 2 =pi -p 2 - s u (20) 
a = (2/t)(v lPl -R + v 2 p 2 -R), (21) 



R = qi + <?2 r = q 2 - q x . (22) 

In particular, the off-shell extension of P x has been chosen in eq. (19) to be more 
symmetric compared to the case considered in ||, and includes two coefficients, V\ 
and v 2 , which are equal to unity at the pion pole. Even with this new extension, 
the polarization vectors still satisfty the requirements 

e «. gi = e (i) -g 2 = 0, 
e «. e (j) = Si, (23) 

for physical polarizations. These relations hold for all positive si and s 2 , whether 
or not they are equal, and for V\ and v 2 chosen as specified in eq. (20). 

The projector n 11 = (n + , n~, n±), using the fact that in the "brick wall" frame 

C W = (0,0, e $?) e i l ) -q ± = 0, (24) 



takes the form 

n 



" ^iPi ± + ^ 2 ± + |(gi ± + g 2 ± ))/v / =i^, 



aq± 
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Notice that only the transversal part of is imaginary. The advantage of using 
such a frame is due to the fact that the spectral densities can be expressed as 
rational functions of the variables Q 2 ,s,Si,s 2 . They are also polynomial in Q 2 
and symmetric in si and s 2 , as expected from the time reversal invariance of the 
scattering process. This provides a drastic simplification of the calculations. It is 
not hard to show that the same simplifications occurs in the case of the pion form 
factor, once it is treated in the same reference frame. 
The final expression for /\( l2 )P ert i s given by 

A (i2 W 5 8Q 2 (s-2Q 2 )R 12 (Q 2 ,s, Sl ,s 2 ) 

3(2vr) 31 (4Q 4 - Sl s 2 f{2Q 2 s - Sl s 2 ) 

8Q 2 (u-2Q 2 )R 12 (Q 2 ,u, Sl ,s 2 ) 



+ 



(4Q 4 - Sl s 2 f(2Q 2 u - Sl s 2 ) 



(26) 



where 
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Ru(Q 2 , s, si, s 2 ) = J2 b n (s, St, s 2 )Q 2n (27) 

n=0 



is polynomial in Q 2 , and with a similar expression for R\ 2 (Q 2 , u, s\, s 2 ) obtained 
replacing s with u in eq. (p7|). This second term comes from the crossed diagram 
of Fig. la, where the two photon lines are interchanged. 

The coefficients b n (s,Si,s 2 ) and b n (Q 2 , u, Si, s 2 ) are symmetric in si and s 2 . 
Their expressions are given by 

bo = -slsl 

b\ = 8s^s 2 (s 1 + s 2 ), 

b 2 = -4sfs2(5ssi + 5ss 2 + 14siS 2 ), 

6 3 = lQsjsl(s 2 si + s 2 s 2 + 14ssis 2 + 2sls 2 + 2s 1 sl), 

b 4 = -16s^2(12s 2 + 15ssi + 15ss 2 + 4siS 2 ), 

65 = 128s 1 s 2 (2s 2 s 1 + 2s 2 s 2 + Ass 1 s 2 + sjs 2 + s 1 sl), 

b 6 = -192sis 2 (4s 2 + ssx + ss 2 + 2s 1 s 2 ) 

67 = 256s(ssi + ss 2 + 2sis 2 ), 

b 8 = 256(-ssi - ss 2 + sis 2 ). 

(28) 

Using the fact that at the double pion pole 

(n- Pl )(n-p 2 ) = s(s + t)/(-t), (29) 
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and inserting this result into eq. (11), we can write down a local duality sum rule 
for Hi + H 2 as 



f 2 (H 1 (s,t) + H 2 (s,t)) = H ds x rds 2 p^ ert ( Sl ,s 2 ,s,t), (30) 

Jo Jo 

with a spectral density given by 



P 



(12)pert 



JLA^ ert ( Sl ,s 2 ,s,t). (31) 



s(s + t) 



By a simple power counting in Q 2 , one can immediately notice that the asymp- 
totic behaviour of the spectral functions, as given by (p6[) show a suppression by 
one additional power of t (t ~ — 2Q 2 , specifically 0(l/t 2 )) compared to the per- 
turbative result predicted by the dimensional counting rules This suppressed 
power-law behaviour, as discussed in [Q, is similar to the one obtained for the pion 
form factor and is not related to any gluonic exchange since the sum rule starts at 
0(a°). The contributions of the power corrections (Figs. 2,3), as we are going to 
discuss next, start at 0(1/Q 6 ) and can be arrested at 0(1/Q 8 ). It turns out that 
a calculation of the power corrections, up to this order, can be done without a full 
evaluation of such diagrams (see appendix B). 



4 Modified Borel transform and gluonic correc- 
tions 

One additional complicacy which appear in the case of 4-point functions, compared 
to 2 and 3 point functions, is related to the fact that the usual definition of the 
Borel transform, so largely employed in the sum rule context, is no more valid. The 
reason is quite elementary and is essentially due to the fact that only a finite region 
of analiticity is allowed for the correlation functions. At fixed angle, such region can 
be depicted as a circle (see Fig. 4) and application of the ope is realistically possible 
only if the bigger radius of the contour is large enough ||. The ope then descibes 
the behaviour of the correlation function exactly on the contour (the far left side 
of it, also termed "the deep euclidean region"). Here we are going to assume that 
only the first lowest dimensional operators are relevant for our needs. Assuming a 
validity of the ope, implemented with only the lowest dimensional condensates, we 
can proceed to the evaluation of such non-perturbative corrections. The coefficients 
(Feynman integrals) of such lowest dimensional operators are, however, analytic 
functions of the pion virtualities inside the circle in Fig. 4 and can be reexpressed 
in the form of contour integrals in the si, s 2 complex planes. Next we are going to 
show that such integrals can be approximated by their double discontinuities parts 
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along the cut from zero to A 2 in the s\ and s 2 planes, integrated on a square of size 
A 2 . To be specific, let's consider in fact the integral operator 



dPl f dp 2 2 /M 2_ p 2 /M 2 



(2vri) 2 Jc Ml Jc Ml 



(32) 



where C denotes any contour that encloses the branch cut in the p\ and p\ planes 
(the shadowed region in Fig. 4), from zero to A 2 , the radius of the region of analitic- 
ity for the correlator of 4-currents, and let's act with it on any scalar integral of 
the form T(p\,p\, s, t), defined by a Cauchy representation of the form 

2 2 1 f f T(S%,S2,S,t) 

«■ «) ~ i <«*/. («,_ A) („ _ pi) ■ < 33 > 

with a possible choice of 7 depicted in Fig. 4. If we let A 4 (si, S2, s, t) denote its 
double discontinuity across the cut, then, using eq. (|32|), we easily obtain 

T 4 (M 2 ,M 2 , S ,t) = B(pj^ Mlpl^ M 2 )T 4 (p 2 ,p 2 , S ,t) 

= £ 2 d 8l £ 2 ds 2 A(s u s 2 , s, t)e- s ^ M ^ s ^ 

(l _ e -(A 2 - Sl )/M^ ^ _ e -(A2- S2 )/Mf ) (34) 



X 



Eq. (|34]) shows that the new transform, acting on each coefficient of the power 
corrections, gives as a result a double integral along the virtuality cut from zero 
to A 2 , times an additional exponential suppression at the edge of the integration 
region. Notice that such additional factors have been introduced into the definition 
of the transform eq. (|32|) in order to ensure that the crossing between the two 
contours C and 7 in Fig. 4 in eqs. (|33|) and (^) does not induce any singularity 
for Si,s 2 = A 2 . We also notice that the influence of such additional suppression 
on the explicit expression of the sum rule is likely to be modest, since the usual 
values of the Borel mass, obtained from the stability analysis of the sum rule for 
the pion form factor, are of the order of 1 Gev 2 , while A 2 is supposed to be higher. 
Therefore, the new definition of the transform compares favourably well with the 
usual one given by |l| 

B(Q 2 - Ml) = Urn Q2 ,_ r^iQ'n-^T- (35) 

In the case A 2 — > 00, B(Q 2 — > Mf ) and ( |35l) coincide. From eq. (0) we come to 
the conclusion that, from each Feynman integral which appears in the evaluation 
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of the power corrections, the modified Borel transform picks-up only the part of it 
integrated along the cut, and therefore we are allowed to neglect all the complex 
contributions in eq. ([33|). The method we implement here is a generalization of the 
one adopted in the simpler case of three-point functions ||. The next point to ad- 
dress in the calculation of the power corrections to the diagram in Fig. 1 is which 
insertion of condensates should be included. Our arguments on this last point are 
entirely based on the similarity between the pion form factor and Compton scat- 
tering. As above, in our kinematic condition, characterized by s,t, u of moderate 
size, all of them in the physical region, and at fixed s/t, operator insertions on the 
quark line connecting the two photons in Fig. 1 are in clear disagreement with this 
similarity picture. By the same token, all the operatorial insertions which lower 
the virtuality of the quark line connecting the two photon lines in the diagrams 
of the power corrections will not be kept into account. A typical example is given 
in Fig. 2b, where the cutting of four quark lines and the specific attachment of 
an internal gluons guarantees continuity of the momentum flow in the loop, but 
lowers the virtuality of the quark line connecting the two photons at the top. Now 
we come to discuss the practical implementation of the method. In the case of the 
power corrections, the general tensor operator G^Gap has been projected onto its 
corresponding scalar average ((G 2 ) = (G a ^ v G a ^ LU ) by the formula 

(0\G^G a p\0) = ^(0\ G ; a G a ^\0}(g, agi/ p - g^ 9va ). (36) 

where l c is the unit 3x3 color matrix. We also assume that 

A Q _ g 2 
Air 

with A a (a = 1, 2, ...8) denoting the standard Gell-Mann matrices. Eq. (^) can be 
also rewritten in the form 

(0|G^G^|0) = ^(OIG^G^IO)^^ - g,pg ua ) (38) 

by using standard properties of the Gell-Mann matrices. 

As usual we work in the Fock-Schwinger gauge defined by the condition 

zM M (z) = (39) 

and we employ a simplified momentum-space method of evaluation of the Wil- 
son coefficients discussed in Ref. |J. In the case of massive quarks, the fermion 
propagator in the external gluonic background is given by |J 

_ j> + m 1. (7" tf P G aP -m la G a ^ p ) 
pz _ m A 2 yp A — m A ) A 

| 7r 2 (G 2 )m £>{m+ fi) 
(p 2 — m 2 ) 4 



G,u = g s G a - a s = ^- (37) 
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which simplifies considerably in the massless limit. Notice that in this limit the 
third term of eq. (|40|), corresponding to the insertion of two external gluon lines 
vanishes and only diagrams with a single gluonic insertion on each fermion line 
(the second term of eq. (Plf)) survive. 

In order to better illustrate our approach (see also the discussion contained in 
Appendix B), let's start considering the scalar integral 

J(pl,pl,s,t,m 1 ,m 2 ) 

d A k 



2 



(pi — k) 2 {k 2 — m 2 ) 2 {{p 2 — k) 2 — m 2 ) 2 (pi — k + qi) 

J(pl,pl,s,t,m 1 ,m 2 ) (41) 



dm 2 dm 2 , 
where we have defined 

~(Pi>P 2 i 55 i) 2) J ^ 2 _ m 2) (Qr> 2 _ fcy _ m 2 ) (j>i — k) 2 (pi — k + qi) 2 ^ ^ 
The dispersive part of eq. ( |42| ) can be approximated as 



J{p 1 ,P2,s,t,mi,m 2 ) = —7-5 / ds x / ds 2 ^- ^ (43) 

47T 2 JO JO (Si -Pi)(s 2 -pi) 



with a spectral density given by 

A scalar ( Sl , s 2 ,s,t,m x ,m 2 ) 



, , ]iJ 6 + (k 2 -m 2 )6 + ((p 2 -k) 2 -m 2 )6 + ((p 1 -k) 2 ) 



(pi-k + q 2 ) 



The dispersive part of the scalar integral associated to eq. fl4T|) when we set mi 
and m 2 to be zero can be written down in the following form 



d 4 k 



J{pl,p 2 2 ,s,t) = J — i t jo 2 



(45) 



4tt 2 

where we have defined 



k) 2 (k 2 ) {{p 2 - k) 2 ) {pi — k-\- qi) 2 

V 2 l d ^L ds \ Sl -pi ){ s 2 -p 2 2 y (46) 



A(si,s 2 ,s,t) = ^ A scalar (s u s 2 , s, t, mi, m 2 )) (47) 



mi,rri2=0 



Evaluation of ( [44]) gives 

25W(mi,m 2 



A scalar (si, s 2 , s, t, mi, m 2 ) = {-2ixif / (46 



12 



with 



W(m 1 ,m 2 ) = -Am\(Q 2 ) 2 + 2m 2 2 Q 2 s + A(Q 2 ) 2 s + 2m 2 1 Q 2 s 1 -2m 2 2 Q 2 s 1 

-2Q 2 ss 1 + 2m\Q 2 s 2 - 2Q 2 ss 2 - m\sis 2 + ssi« 2 (49) 



and with S given as in eq. (|17|). Inserting this result in eq. ([48]) and using fl47|) we 

get 

-327r 4 z(Q 2 ) 2 ( gl - s)(-4(Q 2 ) 2 + 2Q 2 gl - ggi + 2Q 2 s 2 ) 
A(Sl ' S2 ' t} " ,3 (2Q2 _ Sl )3 ( _ 2 g 2 + S2) 3 • (50) 

It is easy to show that the method works perfectly well for values of m x and m 2 
small compared to Q 2 . In fact the spectral functions so obtained are easily shown 
to be regular when the virtualities of the two external pions vary independently in 
the (0, A 2 ) interval. This fact insures are that no poles are generated. However, 
even though the calculation scheme we have developed reduces to few steps, its 
practical implementation requires considerable effort. 

The 3 direct diagrams in Fig. 3 have been evaluated by this method. An 
asymptotic expansion of each coefficient has been made. We get 

A (i2)3a 5a 5 2 (s - 2Q 2 )Q\s - s l ) 2 {2Q 2 s - Sl s 2 ) , * 

91 54tt 2 ^ ; s 3 (2Q 2 - Sl f{2Q 2 - s 2 f 1 ' 

(12)36 20a s 2 (2Q 2 - g ))(Q 4 - AQ 2 s + s 2 + ss, + AQ 2 s 2 ) 

91 ' 27tt 2 ^ ; s(2Q 2 - Sl ) 4 (2Q 2 - s 2 ) 4 1 } 



a (i2)3c 5a s 2 Q 6 (-4Q 4 + 2Q 2 s - Q 2 Sl + as 1 - 3Q 2 s 2 + ss 2 ) 



91 54tt 2V 1 s 3 (2Q 2 - Sl ) 3 (-2g 2 + s 2 ; 



3 



(53) 



where the average on the external gluonic field has been made by using eq. (|36|). 

Crossed diagrams with interchanged photon lines have been treated with the 
same method discussed in this section. However, a much faster way to obtain their 
contributions to the power corrections is by using the s — > u crossing properties of 
the spectral function A^ 12 ). In terms of the variables Q 2 , s and t this is equivalent 
to the substitutions 



Q 2 



Q 2 



u 



Si + s 2 - s - t; 



si + s 2 - 2Q 2 



SlS 2 

2Q 2 ' 



(54) 
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in eqs. (|5l],[52|,[5i|) Combining the contributions from eqs. (j5Tj,52,|53), and adding 
the crossed contributions we can finally represent the gluonic power corrections to 
the spectral function into the form 



A (12) 
gluon 



where 



r^ uon (Q 2 ,s,s 1 ,s 2 ) 



20a s Q s T° luon (Q 2 ,s, Sl ,s 2 ) 
: 27tt 2 s(2Q 2 - Sl )\2Q 2 - s 2 ) 4 

8Q 6 - 24Q 4 s + 14QV - 2s 3 + 7Q i s l 
+2Q 2 ss 1 - 3s 2 s 1 + 17Q 4 s 2 - QQ 2 ss 2 



(55) 



s 2 s 2 



(56) 



This result shows that power corrections for the new sum rules derived in Ref. || 
are calculable in a systematic form in full analogy to the form factor case. The 
power-law fall-off of such contributions is 0(1/Q 6 ), two additional power of Q 
suppressed compared to the perturbative contribution. 



5 Quark power corrections and p( 12 \ 

Quark power corrections can be calculated by the usual methods |J. Typical 
contributions of this type are depicted in Figs. 2. Let's focus our attention on these 
last diagrams, in which the 4 broken lines simbolize a soft insertion of a matrix 
element proportional to ((tpip) 2 ). In the limit of massless quarks the contribution 
proportional to ((?/>?/>)) are zero. As usual ||, the average of the 4-quark operators 
can be obtained from the standard expression 

(0\rj b ^ d s \0) = ^((^) 2 >(rt^ Q ^ 7 - 6™6 bd 6 ai 5 p5 ). (57) 

Combining all possible insertions which satisfy the requirements of leaving the 
top quark line connecting the two external photons virtual and after adding all the 
crossed terms we get 

10TTa s T^ uark (Q 2 ,s 1} s 2 ) 



A^ guark { Sl ,s 2 ,s,t) 



27Q 2 ss?(-2Q 2 + Sl ) 2 s|(-2g 2 + s 2 ) 2 ' 



(5f 



where 

r quark = {-2Q 2 + s )(2Q 2 s-s l s 2 ){s 2 1 + s 2 2 ){16Q 8 - 8Q 6 s - 8QV 
-16Q 6 si + 12Q 4 ssi + 4Q 4 sj - 16Q 6 s 2 + 12Q 4 ss 2 
+16Q 4 Sl s 2 - 2Q 2 ss lS2 - AQ 2 s\s 2 + AQ 4 s 2 2 - AQ 2 Sl s 2 2 + s 2 lS 2 2 ) 



(59) 
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Adding together eqs. (^6| , |55| , |58D the operator product expansion of such spectral 
function is finally written down in the form 



( -t 



s(s + t) / 

1 ) {M 12 ^ ert (s u s 2 ,s,t)0(s -s 1 )0(s -s 2 ) 



\s(s + t) J 
+a s (G 2 ) 



2l 20 Q 8 T3 luon (Q 2 ,s,s u s 2 ) 



27tt 2 s(2Q 2 - Sl ) 4 (2Q 2 - s 2 ) 4 



, 2„ 7n2v 107r Qs r^(g 2 , gl , g2 ) 2 2 

+ 9 (( ^ } > 27Q 2 ssl(-2Q* + Sl ) 2 4(-2Q* + s 2 ) 2 5{Sl ~ Pl) ^ 2 ~ ^ )} ' 



(60) 



where we have included a factor l/{n-p\ n-p 2 ) in the normalization of the spectral 
density and we have used eq. (p9|) As we already mentioned in the previous section, 
our calculations provide us with a spectral density p^ 12 ^(si, s 2 , s,t) which includes 
an exact evaluation of the lowest order contribution (Fig. 1) and is asymptotic in 
l/Q, for the power corrections. Eq. flBTj ), once it is integrated along the (0, A 2 ) cut 
in each of the two virtualities S\ and s 2 , gives a sum rule of the form 



pfpi j h Jo (si -pi){s 2 -pi) 



where we have included all the contributions to the sum rule under a single integral 
sign. In order to further simplify eq. flBTf ) in a consistent way, as we already 
discussed in the previous section, only 3 out of the 8 coefficients given in eq. (|28"D , 
namely bs , 67 and b$ are needed. By doing this approximation in ( |60|) and using 
eq. (126]) it is possible to write down the complete asymptotic spectral function, 
which includes the lowest dimensional power corrections as 

p^ 2 \ Sl ,s 2 ,s,t) 

( ~t \ r 320Q 1 V ert (Q 2 , gl , g2 ) 



5(5 + 1); L 3s(s - 2Q 2 )(4Q 4 - Sl s 2 f(-2Q 2 s + Sl s 2 
20 Q¥'™(Q 2 )S , Sl) s 2 



+« S (G 2 ) 



27tt 2 s(2Q 2 - Sl ) 4 (2Q 2 - s 2 ) 4 



L/ , 7n2N 107ra s r^ fc (g 2 , gl , g2 ) 2 2 

+ (( ^ } ) 27Q^f(=2^ + gl )^(-2Q 2 + , 2 ) 2 ^ (gl ~ Pl) ^ 2 ~ ^ 



(62) 
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where 

T pert = 4 Q4 SSi _ 4 g2 s 2 Si +4 g4 SS2 _ 4 g2 s 2 S2 

-4Q%s 2 - SQ 2 s Sl s 2 + 12s 2 sis 2 (63) 

is the polynomial which appears as a factor in the numerator of the asymptotic 
expansion of the perturbative contribution (the first term in eq. (0)). Notice that 
we still have to act on eq. (|6"ID with the Borel transform eq. fl32|) in order to obtain 
the usual Borel sum rule. In order to show the applicability of the new transform 
and to characterize its action on the power corrections to the spectral function, we 
need some further discussion. 



6 The new Borel sum rule and the pion form 
factor. 

In this last section we discuss two points, namely the possibility to act on the 



sum rule in eq. (|61|) with the modified Borel transform eq. (|32]) in order to get an 
expression which has a similarity with the usual sum rule for the pion form factor 
0, H, at the same time we will argue that the non perturbative power corrections 
are likely to grow much weaker compared to the form factor case in the new sum 
rule for pion Compton scattering. We start from the first issue. 

It is easily observed that the spectral density obtained in eq. ( |5^ ) is regular all 
over the integration region of the two pionic virtualities si ands 2 . This is a nice 
feature of our result, mainly due to the kinematic constraints coming from the fact 
that we are investigating the process at fixed angle and with moderate size s, t and 
u invariants. In fact Q 2 , as given by eq. (|T6| ) is much larger than any of the two 
virtualities even at the edge of the integration region (see Fig. 4) when they are 
approximately equal to (s + t)/A. This fact has an immediate consequence on the 
evaluation of the (modified) Borel transform of the quark power corrections, since 
the only poles which are encountered in the shadowed region in Fig. 4 by the new 
integral transform are those located at p\-,p\ = 0. The evaluation of the residue 
can be obtained, as usual, by expanding the overall factor which multiplies the 
multiple pole l/(PiP 2 ) 3 in the expression of the quark power corrections eq. (|58|) 
up to quartic terms in the virtualities of the two pions. This can be systematically 
done using eqs. (|i~6| , |P71) to reexpress all the powers of Q 2 in terms of si, s 2 and s 
and t. By doing so and by acting afterwards with the modified Borel transform 
fl3~2|) (with Mi = M 2 = M) on the resulting expression we get 

B(pl -> M 2 ,pj -> M 2 )A^ ark (pl,pis,t) 
= na s (F (M, a, t) + F X (M, a, t)e~ x2/M2 + F 2 (M, s, t)e~ 2x2/M2 



(64) 
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where 

F (M,s,t) 



80 (-4M 4 s - 8MV - 2s 3 - 2M 2 st - sH + 2M 2 t 2 + 1st 2 + t 3 ) 

27 MH 2 

8 (2M 2 + s + t)(AM 2 s + 2s 2 - st - t 2 ) 



F±(M, s, t) 



27 MH 2 

(65) 

F 2 {M,s,t) = -™± (66) 

After that this intermediate step is done, we are finally in condition to write down 
a new sum rule for the sum of the two invariant amplitudes H\ and H 2 in the 
helicity base as 



f 2 {^{8,1) + H 2 (s,t)) 



-t 



320Q u TP ert {Q 2 , Sl ,s 2 ) 



X U dSl 'o dS2 3sJs~-2Q 2 )(4Q 4 - Sl s 2 f(-2Q 2 s + s lS2 ) 

+ mT 2as{G) Jo dsi i dS2 S (2Q 2 - Sl )\2Q 2 - S2 y e 

x (l - e -(A a -«i)/^) ^ _ e -(A3- S1 )/M^ 
+ na s ((^) 2 ) (f {M, s, t) + F 1 {M, s, + F 2 (M, s, ^e" 2 ^ 2 ) } 



(67) 



where the expressions of the quark and gluon spectral densities are given in eqs. (|62|). 
Our results has to be compared with its closest one, the sum rule for the pion form 
factor. In this case the usual gluon and quark power corrections give integrals 
which can be explicitely performed after that the standard Borel transform has 
been taken. One gets |]. 0] 

i*SQ P Sty 

f 2 F 7r (Q 2 ) = dsj ds 2 p^{ Sll s 2 ,t)e^ + ^l M2 
Jo Jo 
a s (G 2 } • • n one o ^2 



12ttM 4 

where p n $ has been given in in the form 

pert/ ,\ 
Ptt3 \ s l, s 2,l) 



+ <(^) 2 >f^ra s {l + ^}+ 0(1) +0(« s ) (68) 



3 2 / ( d\ 2 t ( d N " 



? ^ +0 ^ 77- ■ I 7—TTT, (69) 



2tt 4 \\dtj 3 \dt J J ((si + s 2 - t) 2 - 4s 1 s 2 ) 1 /2 
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As we can see from this last equation, the power corrections tend to grow rather 
quickly in the form factor case and invalidate the sum rule at higher values (beyond 
4, 5 Gev 2 ), while the perturbative contributions have a rather fast decay (l/Q 4 ). 
This behaviour of the sum rule sets tight constraints on the range of its validity. 
The behaviour of the perturbative term in (|67D is similar to the form factor case, 
but the power correction are much more stabilized. This fact is related to the new 
angular dependence of the sum rule and to the suppression generated at the photon 
vertex (Fig. 5a). However, in order to give a more precise and definitive answer to 
this last issue, a stability analysis of the new sum rule and a determination of the 
Borel mass M are required. 

7 Conclusions 

We have extended previous work on local duality sum rule for pion Compton scat- 
tering and we have given a complete derivation of the leading spectral function for 
a specific combination of the two invariant amplitudes in the helicity base, H1 + H2, 
showing that pion Compton scattering can be systematically treated as in the pion 
form factor case. The final results for such function turns out to be much more 
complicated compared to the case of the pion form factor, but they still have, at 
moderate momentum transfer, the same power law fall-off, We have seen that the 
sum rule can be organized in the Breit frame of the pion in a rather compact, 
albeit covariant, form. We have also shown that the leading non perturbative cor- 
rections to the sum rule are calculable in a systematic way. For this purpose we 
have found convenient to use some properties of a modified Borel transform which 
acts very simply on all the terms of the sum rule and picks up the residue at the 
pole of the pion virtualities in the quark power corrections. A full discussion of 
the sum rule for the helicity amplitude Hi and the stability analysis of the results 
contained in this work will be considered in a separate work |J . Finally, we expect 
that our results and methods will find applications to a new entire class of elas- 
tic hadron-photon interactions at intermediate energy and momentum transfer. 
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Appendix A. Kinematics. 

Let qi and q 2 be the momenta of the incoming and outgoing photon, respec- 
tively, which are assumed to be on-shell (qf = q 2 = 0). Let also p\ and p 2 be the 
momenta of the incoming and outgoing pion. The external pion states are off-shell 
and are characterized by the invariants s\ = p\ and s 2 —p\- 

We define as usual the Mandelstam invariants 

s = (pi + qi) 2 t=(q 2 -qi) 2 , u = (p 2 - gi) 2 (A.l) 
and the usual relation 

s + t + u = si + s 2 . (A.2) 

It is also convenient to introduce light cone variables for the photon momenta 
as follows 

Qi = qfv + qiV + q 1± , 
v = ^(1,1, ± ), v = -j=(l,-l,0 ± ) 

q 1± . n ± = 0. (A.3) 
In our conventions 

„+ = -L(,,0 + T ,3) tJ= _L(, ; 0_ t; 3 ) _ {AA) 

In the Breit frame of the incoming pion we also easily get 

u = ( P 2-qi) 2 = 2Q 2 -s + ^. (A.5) 

Covariant expressions for qf and q 2 can also be found in the form 

_ (2Q 2 - Sl )(2Q 2 s- Sl s 2 ) 



19 



As discussed in section 3 , a crucial step in the derivation of the sum rule is the 
appropriate choice of the projector n x . Its real part will be denoted as 



n R = (n\n-,-^=) (A.8) 



We easily get the expressions 



(4Q*-4Q 2 s + Sl s 2 ) 
a = ^ (A.9) 



, 2 _ (2Q 2 - s)(2Q 2 - Sl )(2Q 2 - s 2 )(2Q 2 s - Sl s 2 ) 



Note that P 2 , a < 0. 
The relation 



2Q 2 

together with the condition n ■ p\ = n ■ p 2 give 

1 (2Q 2 - Sl )(2Q 2 - s 2 ) 



n-pin-p 2 {2Q 2 - s){2Q 2 s - sis 2 ) 
At the pion pole (si = s 2 = 0, Q 2 = —t/2) we get 

1 _ -t 
n ■ pin ■ p 2 s(s + ty 

which is the result quoted in eq. 



(A.10) 



n-px = (A.ll) 
(s-2Q 2 )(-2Q 2 s + s lS2 ) 



(A.12) 



(A.13) 



Appendix B. Gluonic corrections. 

The contribution to the structure I^ 12 ) for the diagram in Fig. 3b can be 
obtained from the expression 



r^ 2 = (-l/A)G aP G 



pa 



X (2vr) 4 J (p! - k) 2 (k 2 ) 2 ((p 2 - k) 2 ) 2 ( Pl -k + Ql ) 2 

(B.l) 
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whose evaluation proceeds as follows. Following our discussion of section 6 we 
approximate eq. ( |B.1| ) with the double dispersive contribution 

4vr 2 7o Jo (s x -Pi)(s 2 -Pi) 
Introducing the notation 

I[m 1 ,m 2 J(k 2 1 k-p u ...)] 

(B-3) 

d , k f(k 2 , k ■ Pl , ...)6 + (k 2 - ml)5 + (( Pl - k) 2 )6 + ((p 2 - kf - mj) 

(pt-k + qi) 2 

one can give a very compact form to each discontinuity integral. After averaging 
over the gluon field, the spectral function associated to eq. (|B.1|) , in full analogy 
with eq. Q37D is given by 

( d d \ 
A(si, 8 2i s, t) = I ——l[ mu m 2 , f(k 2 , k-pi, ...)] ) (B.5) 



ml,m2=0 



with 



f(k 2 ,k-pi,...) = 64(k ■ n — n ■ pi) 

x (fc 2 /c ■ n — 2k ■ nk ■ pi ~ 2k ■ nk ■ qi 

+k ■ qi"a ■ pi + k ■ ns\ — k ■ np\ ■ qi). (B.6) 

where we have used the relations 

n • q\ = 0, n • n = 0, n • p\ = n ■ p 2 . (B.7) 

Notice that the leading asymptotic behaviour in eq. (p.6|) comes from the term 
containing the factor n-kn-pi, (n-pi) 2 and n-pin-k k-qi which give an asymptotic 
spectral function of the form 

A (I2),asymp 

= ^(G 2 ){I[n -pm-k] (k 2 -2k-pi + si+ I[(n ■ P i) 2 k ■ qi ] 

-3I[n ■ pin - kk ■ qi]}. (B.8) 

Only the leading and first subleading contributions in power of Q, the momentum 
in the "plus" light-cone direction of the incoming pion, are kept. The evaluation 
of the integrals in eq. ( B.8 ) is rather lenghthy. Again, the method to use in this 
case is very similar to the one followed for the scalar case. For completeness we 
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give here the results of the evaluation of 3 main integrals which appear in the case 
of diagram 3a 



I[n ■ pi n ■ k] 
4ttQ 4 (-2Q 2 



s)(-s + s 1 ){-2Q 2 s + sxs 2 ) 



s 3 (2Q 2 - Sl ) 3 (-2Q 2 + s 2 y 



(B.9) 



Y 2 = I[(n ■ Pl ) 2 n ■ k] 

2nQ 4 (-2Q 2 + s)(-s + Sl ) 2 {-2Q 2 s + Sl s 2 ) 

S 3(_2Q2 + Si) 3(_ 2 g2 + S2) 3 > ' 

We also find 

Y 3 = I[n ■ pin ■ k k ■ qi] 

= Y 2 . (B.ll) 

A simple calculation shows that Y\ ~ 1/Q 6 , while Y 2 ,Y 3 m l/Q 4 . Other integrals 
which appear in the power corrections to the spectral function, such as those con- 
taining additional powers of n ■ k are suppressed even more by additional powers 
of Q. Notice also that all integrals we have considered have been evaluated with 
the same method used to derive eq. flSUp. The diagrams in Figs. 3b and 3.c can 
be evaluated similarly. For 3.b we shift with mass parameters the quartic prop- 
agators obtained from the insertion of the external gluon line and proceed as in 
eqs. ( fi6| , |I^ , |50|) . Those terms proportional to additional powers of n ■ k can also be 
evaluated by this method. Their contribution has been included only in the lowest 
order perturbative spectral density eqs. (pp. In the case of the power corrections 
such additional contributions are not going to affect the result to leading and first 
subleading power in 1/Q. 
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Figure Captions 



1. The correlator at lowest order in a s . 
2a, 2b. Typical quark power corrections. 
3a, 3b, 3c. gluonic power corrections. 

The integration contours for the scalar amplitude in the pf plane, showing the 
u = threshold. The shadowed region denotes the area enclosed by the modified 
Borel transform. 

5a. Lowest order perturbative contribution to the spectral density for pion Comp- 
ton scattering. The dashed lines are put on shell. 

5b. Lowest order perturbative contribution to the spectral density of the pion form 
factor. 
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